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Abstract 

We study the late time evolution of flat and negatively curved 
FRW models with a perfect fluid matter source and a scalar field 
having an arbitrary non-negative potential function V ((ft) . We prove 
using a dynamical systems approach four general results for a large 
class of non-negative potentials which show that almost always the 
universe ever expands. In particular, for potentials having a local zero 
minimum, fiat and negatively curved FRW models are ever expanding 
and the energy density asymptotically approaches zero. We investigate 
the conditions under which the scalar field asymptotically approaches 
the minimum of the potential. We discuss the question of whether 
a closed FRW with ordinary matter can avoid recollapse due to the 
presence of a scalar field with a non-negative potential. 



1 Introduction 

Scalar fields currently play a prominent role in the construction of cosmo- 
logical scenarios aiming to describe the structure and evolution of the early 
universe. The standard inflationary idea requires that there be a period of 
slow-roll evolution of a scalar field (the inflaton) during which its poten- 
tial energy drives the universe in a quasi-exponential expansion [1]. Scalar 
fields also arise naturally in alternative theories of gravity which aim to ex- 
tend general relativity, e.g., higher order gravity theories, scalar-tensor and 
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string theories. In higher derivative gravity, due to its conformal relation 
with general relativity [2], scalar fields appear in the Einstein frame with 
a self-interaction, nonlinear potential term which mimics the higher order 
curvature properties of the original (Jordan) frame. Generalizations of the 
original Brans-Dicke theory lead to scalar-tensor theories with scalar field 
self-interactions and dynamical couplings to matter. Further generalizations 
can be achieved by considering multiple scalar fields [3]. Similar results can 
be proved for simple bosonic string theories [4]. 

Most of the studies of scalar-field cosmologies with the dynamical systems 
methods are restricted to Friedmann-Robertson- Walker (FRW) models (see 
for example [5, 6] and references therein). On the other hand, there are also 
important investigations in spatially homogeneous Bianchi cosmologies with 
an exponential potential [7], as well as in models containing both a perfect 
fluid of ordinary matter and a scalar field with an exponential potential, the 
so-called "scaling" cosmologies [8]. 

It is important to investigate the general properties shared by all FRW 
models with a scalar field irrespectively of the particular choice of the poten- 
tial. In this paper we study the time evolution of flat and negatively curved 
FRW models with a perfect fluid and a scalar field with an arbitrary non- 
negative potential function V (0) and prove some general results for a large 
class of such functions. More precisely, we first show that almost always di- 
verges in the past. Furthermore, for potentials having a local zero minimum, 
flat and negatively curved FRW models are ever expanding and the energy 
density asymptotically approaches zero. If in addition the potential has a 
unique global zero minimum, the scalar field almost always asymptotically 
approaches zero. Similarly for exponential potentials it is shown that in an 
expanding universe the energy density p of ordinary matter and the Hubble 
function H asymptotically approach zero, while the scalar field diverges to 
+oo. 

The plan of the paper is as follows. In the next Section we write down the 
field equations assuming an arbitrary non-negative potential, as a constrained 
four-dimensional dynamical system. In Section 3 and prove the above men- 
tioned four propositions under mild conditions on the potential. In Section 4 
we discuss the closed universe recollapse conjecture for a k = 1 FRW universe, 
containing a perfect fluid and a scalar field with a non-negative potential. 
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2 Scalar-field cosmologies 



In General Relativity the evolution of FRW models with a scalar field (ordi- 
nary matter is described by a perfect fluid with energy density p and pressure 
p) are governed by the Friedmann equation, 1 

the Raychaudhuri equation, 

^ = ~1( P + 3P + 2 ^~ 2V )> (2) 
the equation of motion of the scalar field, 

+ 3-0 + V'(0) = 0, (3) 
a 

and the conservation equation, 

p + 3(p + p)- = 0. (4) 

(X 

Here V ((f)) is the potential energy of the scalar field and V' = dV/d<f>. We 
assume an equation of state of the form p — (7 — l)p, with < 7 < 2. 

Formally the FRW model with a scalar field is obtained if we add to the 
matter content of the classical Friedmann universe a perfect fluid with energy 
density e := \4> 2 + V and pressure q := \4> 2 — V . However, this fluid violates 
the strong energy condition, i.e., e + 3q = 2<j> 2 — 2V may be negative. It is 
precisely this violation that leads to inflation in the very early universe [1]. 

From Eqs. (l)-(4) we see that the state (a, a, p, 0, <fij G R 5 of the system 

lies on the hypersurface defined by the constraint (1) and the remaining 
evolution equations can be written as a five- dimensional dynamical system. 
In the case of a flat, k — model, the dimension of the dynamical system 
reduces to four. 

1 We adopt the metric and curvature conventions of [9]. Here, a(t) is the scale factor, 
an overdot denotes differentiation with respect to time t, and units have been chosen so 
that c = 1 = 8ttG. 
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In the remaining of the paper we assume that the potential of the scalar 
field is an arbitrary non-negative (at least C 2 ) function. Gradually we intro- 
duce further conditions on the potential that allow us to analyze the general 
properties of the resulted dynamical system. Setting <fi = y, a/a = H the 
evolution equations (2)-(4) for the flat model become 



subject to the constraint 



4> = y, 

y = -3Hy-V (0) , 
P = ~3lpH, 

ii 1 2 7 
H = ~2 V - 2 P ' 



(5) 



3H 2 = p+-y 2 + V^). (6) 
Therefore, the phase space of the dynamical system (5) is the set 

E := | (0, y, p, H)eR 4 :3H 2 = p + l -y 2 + V (0) j . 

3 Flat and negatively curved FRW with an 
arbitrary non-negative potential 

A theorem due to Foster [10] states that for expanding flat FRW models in 
vacuum with a scalar field having a non-negative potential, almost always 
diverges in the past. The following Proposition is a generalization of this 
result in the case we include a barotropic fluid. 

Proposition 1 Let x G £ and let 0~ (x) be the past orbit o/x under (5). 
Then is almost always unbounded on 0~ (x) . 

Proof. Let x G £ be such that is bounded on 0~ (x) . Using the same 
type of argument as in Theorem 1 of [10], one can show that the a— limit set 
of x lies on the plane p = 0, y = 0. From (5) we see that the only invariant 
sets lying on the plane p = 0, y = 0, are equilibrium points of the form 
(4>,y,p,H) : p = 0,y = and X : V (0i) = 0. Taking into account the 
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constraint (6), we see that the past limit set a (x) consists of equilibrium 
points of the form 0, 0, Hi) , where H\ = V (fa) /3. 

We will show that no such equilibrium is past asymptote to an open 
neighborhood of x. To reduce the dimension of the dynamical system (5) we 
use the constraint (6) to eliminate p. The evolution equations become 

i = y, 

y = -3Hy-V (ft, 

H = + \v (fa . (7) 

For an expanding universe, any equilibrium point of (7) can be written as 

q = ^,0, y/Vo/tj , where fa : V (fa) = and V := V (fa) > 0. It is 

straightforward to see that the eigenvalues of the Jacobian matrix of (7) at 
q are 

- ly /W Q , \ {-y/W Q ± ^(W Q -W» (00)) . (8) 

At least two of them have negative real parts and, therefore the center man- 
ifold theorem implies that there exists a 2 or 3-dimensional stable manifold 
through q (see for example [11]). That means that all trajectories asymp- 
totically approaching q as t — > oo, lie on a 2 or 3-dimensional invariant 
manifold. We conclude that all solution curves asymptotically approaching 
q as t — > — oo, lie on a or 1-dimensional invariant manifold. ■ 

If all of the eigenvalues have negative real parts the assertion follows im- 
mediately from the linearization theorem. In particular, if V" (fa) > 0, 
i.e., fa is a local minimum of the potential, then the equilibrium point 

(^fa, 0, a/Vo/3) is future asymptotically stable. This equilibrium solution 

corresponds to a vacuum de Sitter space with a cosmological constant equal 
to 

Unfortunately the above Proposition cannot be applied to potentials hav- 
ing a zero local minimum V (fa), because for Vo = 0, the eigenvalues (8) have 
zero real parts. However, a large class of potentials used in scalar-field cosmo- 
logical models have a zero local minimum. Examples of potentials belonging 
to this class are polynomial potentials of the form V (0) = \<p 2n , generalized 
and logarithmic potentials studied by Barrow and Parsons [12], or potentials 
in conformally related theories of gravity, for example 

V(fa = Voo(l-eV^) 2 (9) 
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which arises in the conformal frame of the R + aR 2 theory [2, 13]. To study 
the late time evolution of these models, we adopt a different approach. 

In the following discussion we assume that the potential function of the 
scalar field is non-negative with a local minimum V (0) = 0, but otherwise 
arbitrary. For simplicity we present only the flat case, k = 0, but the results 
can be extended to the case of negatively curved FRW models (cf. Remark 
after Proposition 3). 

The equilibrium points of (5) are (0, y, p, H) : p = 0,y = and : 

V (0) = 0. Using the constraint (6) and the assumed V m i n = we see 
that an equilibrium point is the (0, 0, 0, 0) . Physically, the equilibrium point 
(0,0,0,0), corresponds to an empty universe. Linearization of (5) near the 
equilibrium point (0, 0, 0, 0) shows that the eigenvalues of the Jacobian ma- 
trix at that point have zero real parts (the assumptions on the potential 
imply that V (0) = at = and also that V (0) > 0), consequently the 
Hartman-Grobman Theorem does not apply [11]. Therefore, we cannot con- 
clude about the stability of the equilibrium. However, looking at the system 
(5), one arrives at the conclusion that for an ever-expanding universe, H (t) 
and p (t) are decreasing functions and (t) oscillates around the minimum of 
the potential with a decreasing amplitude due to the damping factor 3H(p. 
These intuitive arguments are made more precise in Propositions 2 and 3 
below. 

We first show that an initially expanding flat universe remains ever- 
expanding. In fact, the set {(<f),y,p,H) e M 4 : H = 0} is invariant under 
the flow of (5) with the constraint (6). Therefore, the sign of H is invariant. 
(If the sign of H could change, a solution curve starting with say, a positive 
H, would pass through the point (0, 0, 0, 0) which is an equilibrium solution 
of (5), thus violating the fundamental existence and uniqueness theorem of 
differential equations). We conclude that if H (t ) > 0, then H (t) > for all 
t > to- 
rn the following Proposition the potential function has a local minimum, 

V (0) = 0, but we do not exclude the existence of other local minima with a 
higher value of V. 

Proposition 2 Suppose that V > and V (0) = = 0. Suppose 
also that if A CM is such that V is bounded on A, then V is bounded on A. 
Then lim^ +00 p = = lim^ +00 y. 

Proof. Consider the trajectory passing through an arbitrary point x = 
(0, y, p, H) G £ with H > at t — to- Since H (t) is decreasing and positive, it 
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follows that Hindoo H (t) exists and is a non-negative number 77; also, H (t) < 
H (to), for all t > to- We then deduce from (6) that each of the terms p, \y 2 
and V is bounded by 3H(t ) 2 . Let A = {0 : V (0) < 3H(t ) 2 }. Then the 
trajectory is such that stays inside A. 

From (5) we get 

" f t {^y 2 + \p) dt = H(t)-H (t ) , 
and taking the limit as t — > 00, we obtain 

2 j[ {y 2 + lp)dt = H(t )- V . 

Therefore, 

/•oo 

(y 2 (t)+7p(t))rft<oo. (10) 



/ 



't 

In general, if / is a non-negative function, the convergence of / (t) dt does 
not imply that lim^oo / (t) = 0, unless the derivative of / is bounded. It can 
be shown that lim^oo / (t) = 0, provided that /' is bounded from above. 2 
In our case, 

j t {y 2 + ip) = 2yy + ip = -QHy 2 - 2 y V (0) - 3 7 2 p# < -2yV (0) . 

As we already remarked, y is bounded; also, by our assumption on V, 
V ((f)) is bounded. We conclude that the derivative of the function y 2 + 7p 

2 If / > 0, / < k, where k is a positive constant and J^°° fdt < +00 then, 
lim t ^ +oc / (t) = 0. Here is the proof: Suppose to the contrary that there exists e > such 
that one can find arbitrarily large t such that / (t) > e. Choose such a t and then define 
s < t so that t — so = § • One has 



Vs < i, I f (s) ds<(t-s)k=> f (s) > f (t) -(t-s)k 

J S 



J 

J Si 



* / (s) ds>(t- s ) f (t) ~\k(t- sof > 



The assumption J^°° fdt < +00 implies that lim s ^ +00 f s + °° fds = 0. This contradicts 
the fact that one can find inside any interval (s, +00) a subintcrval (so,t) such that 
Slf(s)ds>^ k . 
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is bounded from above and therefore, (10) implies that lim^ 00 y 2 (t) = 
and lim^oo p (t) — 0. ■ 

The conditions in Proposition 2 are satisfied by a large class of potential 
functions, for example by the generalized and logarithmic potentials studied 
in [12]. In particular, this class includes potentials which do not approach 
zero as — > oo (cf. equation (9)). 

Remark. It is easy to conclude that lim^+oop = 0, without any as- 
sumption of boundedness of V. We only have to consider the case rj > 
because if rj = then (6) entails that lim^oo p = 0. By (5b), p is decreasing, 
thus 3\im t _> +OQ p := p . We shall show that p = 0. Suppose that p > 0. 
Then from (5b) again we get p < — 3?7Po7- Thus 

p(t) - p (t ) = p(s)ds < -3?7po7 (* - *o) , 

J to 

which implies that p becomes negative for large t, a contradiction. 

If we assume in addition that V ((f)) > for > and V ((f)) < for 
(f) < (thus, is the only local minimum of V), then we can conclude that, 
if the initial value of H is less than^/V (oo) /3, then lim t ^ +QO is equal to 0. 
More precisely, we have the following result. 

Proposition 3 Suppose that V ((f)) > for > and V ((f)) < for < 0. 
Then, under the assumptions of Proposition (2), lim t ^ +oo exists and is 
equal to +oo, or — oo or 0. 

Proof. We know that 31im t ^ 00 i/(t) = rj. If rj — 0, then from (6) 
we obtain lim^oo V ((f)) = 0. The assumption on the potential implies that 
lim^oo = 0. So suppose that rj > 0. From (6) we get lim t ^ +00 V (0 (t)) = 
3r) 2 . Thus there exists t' such that V (0) > 3i] 2 /2 for all t > t'. It follows 
that cannot become zero for some t > t' because V (0) = 0. Thus, has a 
constant sign for all t > t'. 

Suppose that > for all t > t'. Since V is increasing as a function of 
on (0, +oo), we have lim^+oo V (0 (£)) = 3rf < lim^ +QO V ((f)). We consider 
two cases: 

(i) If lim t ^ +00 V ((f) (t)) = lim^+oo V (0) then obviously lim t ^ +00 = 
+oo. 

(ii) If lim^+oo V ((f) (t)) < lim^+oo V (0) then there exists > such 
that lim t _^ +00 V ((f) (t)) = V (0). Since V is continuous and strictly increas- 
ing, it follows that lim^+oo = 0. From (5b), taking into account that 
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lim^ +00 y — and H is bounded, we deduce that, lim^+oo y = —V ((f) < 0. 
Hence, there exists t" > t' such that for all t > t", y < -V (0) /2. But this 
entails immediately that 




i.e., y (t) takes arbitrarily large negative values as t increases, which is not 
possible because lim^ +00 y = 0. 

Thus, if > for all t > t' then lim^ +00 = +oo. Similarly, the case 
< for all t > t' leads to \im t _ +OQ = — oo. ■ 

It follows that if initially, 3H 2 (to) < min{lim^ +00 V (0) , lim^_oo V (0)}, 
then lim^ +00 H (t) = 0. Indeed, we know that lim^+oo is equal to or +oo 
or -oo. If, say, lim^ +oo = +oo, then by (6), 3rf = lim t ^ +00 V (0 (t)) = 
lim^+oo V (0) > 3H 2 (t ). This is impossible, since H (t) is decreasing and 
thus H (to) > rj. Likewise, lim t ^ +00 = — oo leads to a contradiction. Hence, 
lim t ^ +00 = and this implies that lim t _^ +00 V (0 (t)) — and, again by (6), 
lim t ^ +00 if (t) = 0. 

We use as an example the potential (9) to visualize the above result. This 
potential has a long and flat plateau. For large values of 0, the potential V is 
almost constant, V^, = lim^+oo V (0) , thus V has the general properties for 
inflation to commence. According to Proposition 3, if the initial value of H 
is less than l^, then 0, H — > as t — > oo. For initial values of if greater than 
V'oo, may diverge to +oo, i.e., the scalar field may reach the flat plateau of 
the potential. 

Remark. The above results can be extended to the case of negatively 
curved FRW models, k = — 1. Setting x — 1/a the evolution equations (2)- (4) 
become 



x 



-Ex 



y 



p 



-3Hy-V'(<f>), 
-3 1P H, 



H 



(11) 



subject to the constraint 




(12) 
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With slight modifications we can repeat the same type of arguments to show 
that Propositions 2 and 3 still hold and, in addition, \im t ^ +OQ x = 0. 

Up to now, we have assumed that the potential is non-negative and has a 
minimum. Therefore, the above three Propositions do not apply directly to 
the important case of an exponential potential, vis. V (0) = Voe~ A ^, with Vo 
and A positive constants. It is already known that the state corresponding 
to p = 0, y = and — > oo is an equilibrium of the projection of the system 
(5) on the — y plane. This result has been established using phase plane 
analysis (cf. [14]) on the reduced two-dimensional dynamical system. Using 
similar arguments as in Propositions 2 and 3, we complete our analysis by 
showing the global result that, for expanding flat models with an exponential 
potential, — > oo as t — > oo. More precisely, we have the following result. 

Proposition 4 Let V be a potential function with the following properties: 3 

1. V is non-negative and lim ( ^_ ( ._ 00 V (0) = +oo. 

2. V is continuous and V (0) < 0. 

3. If A CK is such that V is bounded on A, then V is bounded on A. 
Then lim t ^ +00 y = = lim t ^ +00 p, and lim t _» +00 = +oo. 

Proof. It is easy to see that the set p > is invariant under the flow 
of (5), therefore p is nonzero if initially p (to) is nonzero. Now from (4) it 
follows that H is never zero, thus it cannot change sign. Hence, H is always 
non-negative if H (to) > 0. Further, H is decreasing in view of (5d), thus 
3 lim t ^ +oc H = r] > 0. 

We know from (5d) that 



As in Proposition 2, ^ (y 2 + 7p) < — 2yV ((f)). From this and (13) we deduce 
that lim t _ >+0O y = lim t ^. +0O p = 0, as in Proposition 2. 

We now show that lim^ +00 = +oo as in Proposition 3. It follows from 
(6) that lim t _ i+OQ V (0) = 3i] 2 . Since V is strictly decreasing as a function 
of 0, we have V (0) > \im^ +00 V (0) for all 0, thus \im t ^ +00 V (<p (t)) > 
lim^+oo V ((p). We consider two cases: 

(i) If lim^+oo V ((f) (t)) = lim^ +oo U(0) then obviously lim t ^ +00 = 
+oo. 

3 Note that V ((f)) = Vo e ~ Al * has all these properties. We do not exclude potential 
functions with lim«^ +00 V ((f)) > 0. 




+oo 



(y 2 + 7 p) dt — H (to) -r)< +oo. 



(13) 



to 
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(ii) If lim^+oo V ((f) (t)) > lim^+oo V ((f)), then there exists a unique 
such that lim t ^ +00 V ((f) (t)) — V ((f) . It then follows that lim^ +00 (t) = 0. 
From (5b) we deduce that lim t _ >+0O y = —V' ((f) > 0; thus, there exists 
t' such that y > —V ((f) /2 for all t > t'. From this we get that y (t) — 
y (t f ) > —V (0) (t — t') /2, which is not possible because lim t ^ +00 y — 0. We 
conclude that lim^+oo = +oo. ■ 

If in addition the potential is such that lim^ +00 V (0) = 0, then we 
conclude that H — > as £ — > oo. This is obviously true for an exponential 
potential. 

4 Discussion and the k = +1 FRW model 

In this paper we have investigated the qualitative properties of flat and neg- 
atively curved FRW models containing a barotropic fluid and a scalar field 
with an arbitrary non-negative potential. For potentials having a unique 
minimum V (0) = 0, we have shown that in an expanding universe, the en- 
ergy density p of ordinary matter, the Hubble function H and the scalar field 
asymptotically approach zero. This result was rigorously proved without 
referring to the exact details of the potential. 

However, the situation is more delicate for positively curved FRW mod- 
els with a scalar field. In the present state of the universe the scalar field 
oscillates around the minimum of the potential and it is unobservable, in the 
sense that the energy density p of ordinary matter dominates over the energy 
density e of the scalar field, hence we expect an almost classical Friedman- 
nian evolution. Using (3) it is easy to see that the energy density of the field 
satisfies 

e = -3H(j) 2 

i.e., in an expanding universe e is a decreasing function of time. Since the 
energy density of ordinary matter also decreases, it may happen that in a 
future time, e be comparable to p. In particular, for closed (k — 1) models, 
once the scale factor reaches its maximum value and recollapse commences 
i.e., H < 0, the term 3H(f> in (3) is no longer a damping factor, but acts 
as a driving force which forces the field to oscillate with larger and larger 
amplitude. If this be the case, the repulsive effect of the cosmological term 
may drastically change the evolution of a classical FRW model. 
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These intuitive ideas show that it is of interest to study the late time 
evolution of the dynamical system (11) with k — +1. The main problem we 
are faced with is the following: Can a closed universe filled with ordinary 
matter and a scalar field avoid recollapse? In the following we give a par- 
tial answer to this question for an arbitrary non-negative potential having a 
unique minimum V (0) = 0. 

A closed Friedmann universe is considered almost synonymous to a rec- 
ollapsing universe. This is mainly due to our experience with the dust and 
radiation filled Friedmann models usually treated in textbooks. That this 
picture is misleading follows clearly from an example found by Barrow et al 
[15] according to which an expanding Friedmann model with spatial topology 
S* 3 satisfying the weak, the strong, the dominant energy conditions and the 
generic condition may not recollapse. Thus the problem of recollapse of a 
closed universe to a second singularity is not trivial already in the Friedmann 
case. 

Among the known global results concerning the closed universe recollapse 
conjecture, let us recall that a closed universe recollapses provided that the 
strong energy condition (SEC) is satisfied and there exists a maximal space- 
like hypersurface E, i.e., the expansion is zero on E (cf. [16]). Therefore, in 
FRW models it suffices to check if the scale factor has a maximum. 

In the following we assume an equation of state of the form p = (7 — l)p, 
7 G (|, 2) . 4 Equation (4) can be integrated to give 

p = const x a" 37 . (14) 

Assume that at time t (now) the values of and <fi are very small in the 
sense that 

£ o \i>l + v o < Po and 2$ - 2V < p + 3p (15) 

so that the total stress-energy tensor satisfies the SEC. This is a plausible 
assumption since the scalar field is unobservable in the present universe. We 
write (3) (in first-order approximation) as 

<j> + 3H<p + m 2 = (16) 

where m 2 := V" (0) . This is the equation of motion of an harmonic oscil- 
lator with a variable damping factor 3H. For a slowly varying function H 

4 The range of 7 is chosen in accordance to the conditions for recollapse of Barrow et 
al [15]. 
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this equation can be solved using the Kryloff-Bogoliuboff [17] approxima- 
tion. We find that the amplitude of the scalar field varies as ~ a -3 / 2 . Since 
the amplitude of has the same time dependence and in our approximation 
e = \4> 2 + |m 2 2 , it follows that 

e ~ a" 3 . (17) 

Comparing this with the time dependence of the density p ~ or 37 we arrive 
at the following picture for the evolution of these universes. If 7 < 1 the 
initial conditions (15) imply that 

e = -<p 2 + V < p and 20 2 - 2V < p + 3p, 

for all t >t , that is, the SEC on the total stress-energy tensor is satisfied for 
all t > t . Hence the universe follows the classical Friedmannian evolution 
and has a time of maximum expansion. Therefore, with initial conditions 
(15) the approximation (16) suggests that if 7 < 1, the universe recollapses. 
The case 7 = 1 is particularly interesting since it corresponds to a dust- 
filled universe which perfectly approximates our Universe. For 7 > 1, since 
p decreases faster than the energy density e of the scalar field, e eventually 
dominates over p. Hence this model may, or may not have a time of maximum 
expansion. 

The above analysis may serve as an indication of the richness of possible 
scenarios for expanding closed FRW models with a scalar field. However, 
spurious conclusions may be drawn from an approximate analysis, especially 
if it is referred to the asymptotic behavior of the solutions. Similarly, the 
determination of the stability of equilibrium solutions does not reveal the 
complete global behavior of a dynamical system of dimension greater than 
two. A rigorous proof of the closed universe recollapse conjecture may come 
from the investigation of the global structure of the solutions of (11) with 
k = +1. 

Acknowledgments I thank N. Hadjisavvas and S. Cotsakis for fruitful 
discussions during the preparation of this work. 
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